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Abstract 

Network calculus is a powerful methodology of characterizing queueing processes and has 
wide applications, but few works on applying it to 802.1 1 by far. In this paper, we take one 
of the first steps to analyze the backlog bounds of an 802. 1 1 wireless LAN using stochastic 
network calculus. In particular, we want to address its effectiveness on bounding backlogs. 
We model a wireless node as a single server with impairment service based on two best- 
known models in stochastic network calculus: Jiang's and Ciucu's. Interestingly, we find 
that the two models can derive equivalent stochastic service curves and backlog bounds in 
our studied case. We prove that the network-calculus bounds imply stable backlogs as long 
as the average rate of traffic arrival is less than that of service, indicating the theoretical ef- 
fectiveness of stochastic network calculus in bounding backlogs. From A. Kumar's 802.1 1 
model, we derive the concrete stochastic service curve of an 802. 1 1 node and its backlog 
bounds. We compare the derived bounds with ns-2 simulations and find that the former 
are very loose and we discuss the reasons. And we show that the martingale and indepen- 
dent case analysis techniques can improve the bounds significantly. Our work offers a good 
reference to applying stochastic network calculus to practical scenarios. 
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1 Introduction 



Network calculus provides an elegant way to characterize traffic and service pro- 
cesses of network and communication systems. Unlike traditional queueing the- 
ory in which one has to make strong assumptions on arrival or service processes 
(e.g., Poission arrival process, exponential service distribution, etc) so as to derive 
closed-form solutions in queueing networks [JJ, network calculus allows general 
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arrival and service processes. Instead of getting exact solutions, one derives net- 
work backlog and delay bounds by network calculus. Deterministic network calcu- 
lus is mature in theory [2|l iSl B [|3. However, most traffic and service processes 
are stochastic and deterministic network calculus is often not applicable to them. 
Therefore, stochastic network calculus was proposed to deal with stochastic arrival 
and service processes [|7]|- ifTSll ll25l . 

Numerous applications of it have been found in communication networks and even 
in management science, and we cite some of them |fT9l - [|27l . However, few works 
have been made on applying it to multiple access communication networks such as 
802. 1 1 Wireless LANs |[26ll [|25l . In the paper, we take one of the first steps to apply 
stochastic network calculus to an 802.11 wireless LAN (WLAN). In particular, 
we want to address the effectiveness of stochastic network calculus on bounding 
backlogs in 802.11, with the following sub-problems: 

• Under what condition can we derive stable backlogs using network calculus? 

• How to derive the concrete stochastic service curve an 802.1 1 node? 

• Are the derived backlog bounds tight compared with ns-2 simulations? And how 
to improve them? 

We model a wireless node as a single server with impairment service based on two 
best-known models in stochastic network calculus: Jiang's [[T6ll and Ciucu's [fTTl . 
And we make the following new contributions on this topic: 

• We compare Jiang's and Ciucu's model and find that they can derive equivalent 
stochastic service curves and backlog bounds in our studied case. 

• We prove that the network-calculus backlog bounds imply stable backlog as long 
as the average rate of traffic arrival is less than that of service, indicating that 
stochastic network calculus is effective in bounding backlogs theoretically. 

• From A. Kumar's 802.1 1 model, we derive the concrete stochastic service curve 
of an 802.11 node [30] and give the numerical computation methods. From the 
service curve we then derive backlog bounds. 

• We observe the derived bounds are loose when compared with ns-2 simulations. 
However, the martingale and independent case analysis techniques can improve 
the bounds significantly. 

Note that when we prove a statement in this paper, we call it Propositions to differ- 
entiate the existing theorems in the literature (see Proposition [UlSl). 

This paper is organized as follows. In Section [2l we give a brief overview of 
stochastic network calculus. In particular, we present the classic models of Jiang's 
and Ciucu's and we also discuss the martingale and independent case analysis tech- 
niques. In Section [3l we present the network calculus model of a wireless node 
based on Jiang's and Ciucu's model. We compare the two models and find that they 
are equivalent in deriving stochastic service curves and backlog bounds in our stud- 
ied case. We also prove the stability condition by the theory of stochastic network 
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calculus in this section. In Section HI we derive the backlog bounds of an 802.11 
node and the critical part is to derive its concrete stochastic service curve. In Sec- 
tion [51 we compare the derived backlog bounds with ns-2 simulation results under 
Poisson traffic arrivals. In particular, we show that the martingale and independent 
case analysis techniques can improve the bounds significantly. In Section [6l we 
give related works and highlight our contributions. Finally, Section|7] concludes the 
paper and points out some future works. 



2 Stochastic Network Calculus 

In this section, we first review basic terms of network calculus and then cite some 
results of the stochastic network calculus theory used in our paper. Jiang classified 
stochastic arrival curves as the types of ta (traffic amount centric), vb (virtual back- 
log centric) and mb (max virtual backlog centric), and classified stochastic service 
curve as ws (weak stochastic) and sc (stochastic). In this paper, we adopt ta and mb 
arrival curves and the ws service curve, as currently they provides tightest back- 
log boundfl Note that we just say "stochastic service curve" in our paper which 
means the ws one. 



2. 1 Basic Terms of Network Calculus 

We consider a discrete time system where time is slotted (t = 0,1,2,...).A process 
is a function of time t. By default, we use A{t) to denote the arrival process to a 
network element with A(0) = 0. A{t) is the total amount of traffic arrived to this 
network element up to time t. We use A*{t) to denote the departure process of the 
network element with A*{Q) = 0. is the total amount of traffic departed from 
the network element up to time t. Let (^) represents the set of non-negative 
wide-sense increasing (decreasing) functions. Clearly, A{t) e T and A*{t) e J-". 
For any process, say A(t), we define A{s, t) = A{t) — v4(s), for s < t. We define 
the backlog of the network element at time t by 

Bit) = Ait)-A*it), (1) 
and the delay of the network element at t by 

^ As recently known by the network calculus community the bounding probability of a mb 
arrival curve is 1 for linear arrival curve functions, making its usage restrictive. So does the 
sc service curve as it is often derived from an impairment process with the mb arrival curve. 
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D{t)=M{d: A{t) <A*{t + d)}. (2) 
Fig. □ illustrates an example of A{t) and A*{t) with B{t) and D{t) att= 10. 
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Fig. 1. Illustration of A{t), A*{t), B{t) and D{t) 



In deterministic network calculus, A{t) can be upper-bounded by an arrival curve. 
That is, for all < s < t, we have 

A{s,t) < a{t - s), 

where a{t) is called the arrival curve of A{t). 

We say, busy period is a time period during which the backlog in the network 
element is always nonzero. For any busy period (to, t], suppose we have 

A\t)-A*{U)>l3{t-to). 

if the network element provides a guaranteed service lower-bounded by — to) 
during the busy period. We can let to be the beginning of the busy period, that is, 
the backlog at to is zero or A* (to) = Aito). Therefore, 

A*(t)-A(to) >/3(t-to). 



The above equation infers A*{t) > info<s<i [A{s) + /3(t — s)], which can be writ- 
ten as 



A*{t)>A®p{t), (3) 

where ® is called the operator of min-plus convolution and /3(t) is called the service 
curve of the network element. 
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2.2 Stochastic Network Calculus Theory 

We consider a server S (i.e. the network element) fed with a flow A. In practice, 
A's traffic and S"s service are often stochastic, which can not be hard bounded 
by some curves. That is, they can violate the curves but with certain probabilities 
(we call it bounding function here). The theory of stochastic network calculus can 
get probabilistic bounds for backlogs and delays of the server, suppose we can 
characterize A by a stochastic arrival curve and S' by a stochastic service curve. 

In this section, we just consider the derivation of backlog bounds as delay bounds 
are quite similar to the former. We first give some definitions. Then we cite some 
results in Jiang's and Ciucu's models [fT6ll [fTTl and the construction . Lastly, we 
make a brief discussion on them. 



2.2.1 Definitions 

Definition 1 (ta stochastic arrival curve) A flow is said to have a ta (traffic-amount- 
centric) stochastic arrival curve a E with bounding function f E denoted by 
A ^ta< f, OL >, if for all s,t > 0(s < t) and all x > 0, there holds 

P{A{s, t) - a{t -s)> x} < f{x). (4) 

Definition 2 (vb stochastic arrival curve) A flow is said to have a vb (virtual- 
backlog-centric) stochastic arrival curve a E J-' with bounding function f E T, 
denoted by A ~^b< /, a >, if for allt > and all x > 0, there holds 

P{ sup [A{s, t) - a{t - s)] > x} < f{x). (5) 

0<s<t 

We can see that A ~t,fe< f,a> implies A ~ja< f, « since P{A{s, t) — a{t — 
s)> x} < P{supo<,<j[A(s,t) - a{t - s)] > x}. 

Definition 3 (Stochastic Service Curve) A server S is said to provide a (weak) 
stochastic service curve (5 E with bounding function g E T, denoted by S < 
g, /3 > (or just S ~< g, 13 >), if for allt > and all a; > 0, there holds 

P{A^f3{t)-A*{t)>x}<g{x). (6) 

Definition 4 (Leftover Service) Consider a server S provides the ideal service 
curve Pit) with the impairment process I to a flow. Then, during any backlogged 
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period (s, t], the output flow A*{s, t)from the server satisfles 



A*{s,t)>P{t-s)-I{s,t). (7) 
/3(t) — I{t) is the leftover service received by the given flow. 

The definition of leftover service (also called stochastic strict server in [16] ) can 
be applied to many scenarios such as cross traffic and wireless channels. 

Definition 5 (6'-MER) A process A 's minimum envelope rate with respect to 9 (9- 
MER), denoted by p*{9), is deflned as follows: 



p* {9) = Um sup ^ sup log £e^^('''+*) . (8) 

t^oo 9t s>0 

We say that A has an envelope rate with respect to 9 (9-ER), denoted by p{9), if 
p{9) > p*{9). 

Definition 6 ((cr(0), p(^^))-upper constrained) A process A is said to be {a{9),p{9))- 
upper constrained /or some 9 > 0, if for allO < s <t, we have 



I log ^e^^(^'*) < p{9) (t - s) + a{9) . (9) 
9 



We can derive stochastic arrival and service curves from the {a (9), p{9)) -upper 
constrained characterization (Section l23l 

Definition 7 (Average Rate) The average rate of a process A, denoted by a a, is 
deflned as 



ga = iim sup . (10) 

s>0 t 

Definition 8 (Stable Backlog Bound) The backlog B{t) is stable, if for all t, 



EB{t)<C <oo, (11) 

where C is aflnite constant value. We say that the backlog bounds are stable if they 
can derive stable backlogs. 
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2.2.2 Jiang's Model 



Jiang's model deals with vh arrival curves and stochastic service curves. We have 
the following theorems for leftover service curves and backlog bounds. 

Theorem 1 (Jiang's Leftover Stociiastic Service Curve) Suppose a server S pro- 
viding the ideal service curve (3{t) with the impairment process I. If I has a vb 
stochastic arrival curve, i.e., I ^vb< 9,1 >, then the server provides the flow the 
leftover stochastic service curve S ^< g, P > and 



m-m-i{t). (12) 

Theorem 2 (Jiang's Baclilog Bounds) If the flow A has a vb stochastic arrival 
curve A ~^(,< f,Q:> and the server S provides a stochastic service curve 5 ~< 
g, P > to the flow, then the backlog B{t) of the flow in the server at time t satisfles: 



P{B{t) > X + sup [a{s) - p{s)]} < / ® g{x), (13) 
for allt > and all x >0. 



2.2.3 Ciucu's Model 

Ciucu's model deals with ta arrival curves and stochastic service curves. 

In fact, we can derive vb arrival curves from ta arrival curves by introducing the 
function 6{t) = 6 ■ t (5 is an adjustable constant). The following lemma states this. 

Lemma 1 (ta to vb Arrival Curves) Suppose A is a ta stochastic arrival curve, A ^ta< 
f,a >, then A ^vb< f,Ci6> with as{t) = a{t) + 5t and its bounding function 
f{x, S) — Z^^o /(^ + (suppose the sum isflnite). 

The derivations are as follows. 
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P{ sup [A{s,t) — as(t — s)] > x} 

0<s<t 
t 

<J2P{A{s,t)-as{t-s)>x} 

s=0 

t 

= P{A{s, t) - a{t - s) > X + 5{t -s)} 

s=0 

t oo 

< ^ /(x + Sit -s))<Y: fix + kS). (14) 

s=0 k=0 

We have the following theorems for leftover service curves and backlog bounds in 
Ciucu's model. Actually, we can derive these results by first converting ta arrival 
curves to vb ones and then applying Jiang's theorems. 

Theorem 3 (Ciucu's Leftover Stochastic Service Curve) Suppose a server S pro- 
viding the ideal service curve Pit) with the impairment process I. If I has a ta 
stochastic arrival curve, i.e., I ^ta< 9,1 >, then the server provides the flow the 
leftover stochastic service curve S ~< g, (3 > and 



m = m-isit), (15) 

where Jsit) = 7(t) + 6t and gix, 6) = J^'kLo di^ + k5) by definition. 

Theorem 4 (Ciucu's Backlog Bounds) If the fiow A has a ta stochastic arrival 
curve A ^ta< f)Oi> and the server S provides a stochastic service curve S ~< 
g, P > to the fiow, then the backlog Bit) of the fiow in the server satisfies: for all 
t > and all x > 0, 



P{Bit) > a: + sup [asis) - Pis)]} < / ® gix), (16) 

where asit) = ait) + 6t and fix, 6) = J^kLo f{x + k6) by definition. 

Note that Ciucu's can deal with ta arrival curves while Jiang's can not, by introduc- 
ing a 5 > to trades smaller service for finite bounding functions. 



2.3 Computation of Stochastic Arrival/Service Curves 



We will show in this subsection how to calculate stochastic arrival and service 
curves from the iaiO), pi9)) -upper constrained characterization [7J. 
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Theorem 5 (Arrival Curves of {a{6), p{6))-upper constrained) Suppose A{t) is 
{a{9), p{6))-upper constrained, then it has a ta stochastic arrival curve A ^ta< 
f, a >, where 



a{t)=r ■ t 

/(a;)=e^'^(^) -e-^", (17) 

for any r > p{9) and x > 0. And A has a vb stochastic arrival curve A ~^5< 
/, a >, where 

a(t) = r ■ t 

eu(e) 

/(^) = l_e.(p(.)-.) -^"" (18) 
for any r > p{6) and x > 0. 

Note that we have r > p{0) in ta and r > p{6) in vb. And Eq.(fT8l) applies Boole's 
inequality to the bound functions f{x) which are loose in general. 

How to derive stochastic service curves? If we can model the server S with the 
ideal service curve (3 with the impairment process I{t), we can first characterize 
J(t) by vb (ta) arrival curves, and then we use Theorem [T] (Theorem (3]) to get its 
stochastic service curves. 

The following theorem states the relation between 6'-ER and ((t(6'), p(^))-upper 
constrained. We will use it in proving the stability condition of backlog bounds in 
Section [331 

Theorem 6 (6'-ER vs (cr (6'), p(6')) -upper constrained) If the process A{t) has a 
6-envelop rate (9-ER) p{6) < oo, then for every e > there exists cr^(9) < oo so 
that A is {(T^{9), p{9) + e)-upper constrained. 



2.4 Improvement on Bounds 



There are two ways of improving bounds in current literature. One way is to apply 
independent case analysis. The other way is to improve the bounding functions of 
stochastic service curves for time-independent arrivals. 

The first way says that suppose the impairment process of the server S is indepen- 
dent from the traffic arrival process, we can derive tighter backlog bounds using 
independent probability analysis. 
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Theorem 7 (Backlog Bounds under Independent Cases) Suppose the server S 
provides the flow (satisfying A ~^f,< f,a>) the ideal service curve /3{t) with the 
impairment process I ^vb< 9,1 > (thus S ~< g, (3 > where f3(t) = f3(t) — 7(t)j. 
Suppose A and / are independent, we have 

P{B{t) > sup[a(s) - I3{s)] + x} 

<Y^m)-9{k-l))f{x-k) (19) 

fc=0 

where f(x) = 1 — f\x), g{x) = 1 — g{x), and we set 1) = 0. 

This theorem of independent case analysis can be applied to Ciucu's model. How- 
ever, we first need to convert the ta arrival curves to the vb ones with new bounding 
functions /(a;, 5) and g{x, 5) by Lemma [TJ Then we apply the above theorem by 
plugging in / and 5-. 

Another way of tightening backlog bounds is to derive tighter bounding functions 
of stochastic arrival and service curves. Ciucu first proposed to use martingale to 
tighten the bounds for M/M/1 and M / D /I queues [18]. In the following propo- 
sition, we provide a more general result following his idea. The proof is given in 
Appendix- A. 

Proposition 1 (vb Arrival Curves of Time-Independent Process) Suppose A{t) is 
{(J (6), p{6)) -upper constrained. On condition that a{t) = A{t) — A{t — 1) is inde- 
pendent of each t, it has a vb stochastic arrival curve A ^vb< f, ol >, where 



a{t) 
fix) 



-r ■ t 



(20) 



for any r > p{9) + a{9) and x > 0. 



2.5 Discussion on Jiang's and Ciucu's Models 



The key difference between Jiang's and Ciucu's models is: Jiang use vb traffic 
arrival curves while Ciucu uses ta ones variant stochastic service curves. Which 
one can derive tighter backlog bounds? 

In general, ta arrival curves provide tighter bounding functions than vb. Actually, 
^ ^vb< f,0!> implies A ^ta< f-iOL> and the inverse is not true generally. In 
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particular, the bounding function of ta is tighter than that of vb (especially when r is 
close to p{9)) in Theorem[5l But one can not conclude that Ciucu's model is always 
better than Jiang's, as it has looser bounding functions for leftover service curves 
and backlog bounds (see Theorem [3] and Theorem IH). The situation becomes even 
more uncertain when consider time-independent processes and independent A and 
/. Interestingly, we find that two models can derive equivalent stochastic service 
curves and backlog bounds in our studied case (Section [33)) . 



3 A Wireless Node's Network Calculus Model 

In this section, we model a general wireless node by stochastic network calculus. In 
general, we can define one time slot (t = 1) to be any small duration and measure 
traffic amount in any unit (e.g. bits, bytes or packets). 

We consider a wireless node. Let A{t) denote the traffic arrived at the node from 
the application layer. We assume A is (cr^(6'i), p^(6'i)) -upper constrained, which is 
a right assumption for many cases. 

We model the service of a wireless node as an ideal server curve with an impair- 
ment process. In fact, the channel is shared by the other node in a WLAN and 
transmission errors occur due to path loss , fading and coUisons, which contribute 
to the impairment process. Let the channel capacity be c traffic units per slot. The 
departure process A*(s, t) = (3{s, t) — I(s, t) during any backlogged period [s, t], 
where f3{t) = c ■ t is the ideal service curve and / is the impairment process. Since 
J(s,t) < c- (t — s), there exist (T/(6'2) andp/(6'2) so that lis (cr7-(6'2), p/(^^2) ) -upper 
constrained. 

In here, Oi and O2 are adjustable parameters. We will show Section|4]how to calcu- 
late pa(^i), (yA{Oi), Pi (02) and ai{92) for an 802.11 node. 



3. 1 Jiang 's Backlog Bounds 

Because A is (ctaI^^i), Pa (^^i)) -upper constrained, by Theorem [51 A ~^6< f,a > 
where 

a(t) = rA ■ t 
for any > pa{Oi). 
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In the same way, I ^vb< 9,1 > where 



l{t) = ri-t 




1 _ ee2(piie2)-rj) 



(22) 



for any rj > pi{92)- 

By Theorem [H the node provides a stochastic service curve S g, (3 >, where 



for any c> rj. 

Finally, by Theorem[2l we must let a{t) < (3{t), i.e., < c — r/, in order to get 
meaningful backlog bounds which are P{B{t) > x} < f ® g{x). 

We note that f{x) {g{x)) is the decreasing function of (r/). Considering the 
above conditions, we get the following optimal backlog bounds, 

P{B{t) > x} < min [/ ® g{x)] 

9i,02,rA,ri 

subject to 

ta > PA{Oi),ri > pi{02) 

TA + rj = c 



In here, Pa(^i) (P/(^2)) is the function of 9i {62). 



3.2 Ciucu 's Backlog Bounds 

Because A is ((Tyi(6'i), pA (6*1)) -upper constrained, by Theorem [51 A ^ta< f,ct >, 
where 

a(t) = rA ■ t 



P{t) = ic-ri)-t, 



(23) 



^1,^2 >o. 



(24) 



/(x)=e' 



(25) 



for any > pa{Oi). 



In the same way, / ~ja< 9,7 >, where 
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-f{t) = rj-t 

= e^2"^(^2) • e-^2^ (26) 
for any r/ > p/(6'2). 

By Lemma [H we have A ^vb< /, tt<5i > where 



«<5i(^) = i^A + -t 

/(^'^i) = T3-Mr-^"'"' (27) 



for any > Pa(^i) and (5i > 0. In here, we can get the close form of f{x, 5) 
Efclo fi^ + for the particular f{x) in Eq.dlSl). 

In the same way, / ~t,6< g, > where 



l5At) = {ri + 62)-t 
for any r/ > piiQ-i) and ^2 > 0. 

By Theorem [31 the node provides a stochastic service curve ~< ^, /3 >, where 

/3_,,(t) = (c-r,-52)-t, (29) 
for any c > r/ + (52- 

Finally, by TheoremlH we must have 0:5^ (t) < /3_52(t), i.e., r^ + 5i < c — r/ — 52in 
order to get meaningful backlog bounds which are P{B(t) > a;} < / ® gix). We 
note that /(x) ((7(x)) is the decreasing function of 61 (52). Considering the above 
conditions, we get the following optimal backlog bounds. 



P{B{t)>x}< min [f(^g{x)] 

8i,62,Si,02,rA,ri 

subject to 

rA > pA{di),ri > pi{92) 
rA + rj + 61 + 62 = c 

61,62 > 0,61,62 >0. (30) 



In here, Pa{6i) (pi{62)) is the function of 61 (62). 
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3.3 Equivalent Bounds in Two Models 



We find that the two models actually can derive the same stochastic service curves 
and backlog bounds. To understand this, we note that the key difference is the traffic 
model. The following proposition shows that we can derive the same vb arrival 
curves from the two models. 

Proposition 2 (ta vs vb Arrival Curves) If A is a {a{9), p{6))-upper constrained 
process, its vb arrival curve immediately generated by applying Theorem\5^s Eq. dTSj) 
and the one generated by applying Theorem \5\'s Eq. n7\) and then Lemma [7] are 
equivalent. 

Proof: Following the discussions above, the vb arrival curves generated immedi- 
ately by applying Theorem's Eq.(fT8]) are A ~j,fe< /, a > where 

a(t) = r ■ t 

eaie) 

for any r > p{6). 

The vb arrival curves by applying Theorem [5js Eq.(fT71) and then Lemma [T] (con- 
verted by the ta arrival curves) are A ^yb< f,c(s> where 

as{t) = {r + 5)-t 
ea{e) 

/(x,5) = ^— 3^•e-^ (32) 
for any r > p{9) and 6 > 0. 

For the same value of {r+5) in Eq.(l32l). we should maximize S to get tighter f{x,6); 
in other words, we should minimize r and let it to be p{9) . In this optimized case we 
find that Eq.(l3TI) and Eq. (|32l) are in the same form. This establishes the equivalence 
between them. ■ 

This result can be applied to the vb arrival curves of impairment processes. Since 
Ciucu's model can be derived from Jiang's for the single-server case (Section r2.2.3l) . 
the two models can derive the same backlog bounds in the following proposition. 

Proposition 3 (Bounds Equivalence in Two Models) Consider a single server S with 
an ideal service curve (3 and an impairment process I. Suppose the traffic arrival 
process A and the impairment process I are {cr{6), p{6))-upper constrained for 
some 9 respectively, then the vb arrival curves, the stochastic service curve and 
backlog bounds derived by Jiang 's and Ciucu 's model are equivalent. 
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We omit the proof which is got immediately from Proposition [21 

Note: By equivalence, we do not mean it is general for all situations. Actually, 
Ciucu's model can be extended to multiple concatenated nodes while Jiang's can 
not, and they are different models. Even for the single-node case, we only prove the 
equivalence property for linear arrival curves in our studied case. And it is still an 
open problem for more general cases. 



3.4 Stability Condition 

One fundamental question we need to address is under what condition we can de- 
rive stable backlog bounds (i.e., Ei?(t) < oo) by stochastic network calculus. The 
following proposition shows the stability condition. 

Proposition 4 (Stability Condition) Suppose there exist 9-MERs ( 9-Minimum En- 
velop Rates) for the traffic arrival process A and the impairment process I of the 
wireless node for < 9 < 9 where 9 is some constant value, then stochastic net- 
work calculus can derive stable backlogs if 

qa < c — aj, (33) 

where c is the transmission rate of the ideal channel, a a and aj are the average 
rate of A and I defined in Definition^ respectively. 

Proof: 

The proof consists of two phases. First, we show that a a < c — aj can lead to 
ta < c — ri. Next, we show that if < c — r/ then stochastic network calculus 
can derive EB(t) which is less than a finite value. 

We adopt Jiang's model in S ection [311 where A is the traffic arrival process and / is 
the impairment process of the server S, since the two models are equivalent in our 
studied case (see Proposition [3]). We have shown that P{B{t) > x} < f ® g{x) if 
< c — r/ holds. 

From Eq. ((2T[) and ([22]) . we let ei = — Pa{9i) and 62 = r/ — pi{92) for 6'i, 6*2 > 
and ei, e2 > 0. To simplify the arguments, we let 9i = 92 = 9 and ei = €2 = e. 

Thus, < c — r/ holds if 

PAi0)<c~pii9)-2e. (34) 
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From Theorem [6l we can construct the {(j{9), p(^^))-upper constrained characteri- 
zation by letting Pa(^) = Pa(^) + ^ ^^'^ Pii^) = P*ii^) + ^ e > 0, where 
P*a{Oi) and p}{62) are 0-MERs of A and /, respectively. And Eq.dlll) holds if 



PAW<c-pK^)-4e. (35) 
Because p*a{0) exists, applying Taylor's expansion, 



Pa{9) = lim sup ^ suplogEe^^(''*+*) 

t^oo ut s>0 

1 



lim sup — sup logE(l + s + t) + 0(6^ A{s, s + ty)) 
lim sup ^ sup log (1 + eEA{s, s + t) + 0(^M(s, s + tf)) 
lim sup ^ sup [^EA(s, s + t) + 0{e'^A{s, s + tf)]. 



Let e go to 0, 



EA(g,g + t) 

limp4(6') = lim sup = O/i. (3d) 



Similarly, 



limpK^) = a/. (37) 



Therefore, there exists some 9 < 9 so that p*a{9) < a^i + e and p}{9) < aj + e. So 
Eq. dMl) holds if 



< c — a/ — 6e. (38) 
Since e can be arbitrarily small, Eq.(l38]) holds if 

ttA < c — aj. (39) 
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Following the above derivations backwards, we prove that < c — aj leads to 

< C — Tj. 

Next, we prove that stochastic network calculus can derive EB{t) which is less than 
a finite value if r a < c — rj. 

Since f{x) and g(x) are exponentially decreasing functions according to Eq. (|2TI) 
and Eq. (l22l) . we can show that EB(t) is upper-bounded by some finite constant 
value as follows. Note that B{t) is a discrete value in practice (e.g., in bits or pack- 
ets). 



EB{t) = J2 P{B{t) = k + l}-{k + l) 

oo 

<J2P{Bit)>k}.ik + l) 

k=0 
oo 

<Y.f^9{k)-{k + l) 

k=0 

°° k k 

<E(/(bJ) + ^?(M))-(^ + i)<oo. (40) 

fc=0 ^ ^ 



Remarks: Since the proof is based on the theory of stochastic network calculus, it 
indicates that we can get stable backlog bounds by stochastic network calculus on 
the condition that the average arrival rate is less than the average service rate. As 
this condition is very general, stochastic network calculus is effective in theory. 



4 An 802.11 Node's Network Calculus Model 

In this section, we derive the backlog bounds for an 802. 1 1 node. And the key part 
is to derive its stochastic service curve. We use Jiang's model in this section since 
the two models are equivalent in our studied case (Proposition [3]). 

For simplicity, we assume n identical nodes send packets to an AP (access point) 
and they share the wireless channel. All nodes operate in Distributed Coordination 
Function (DCF) mode with RTS/CTS turned off [28]. We assume that transmission 
errors only happen due to packet collisions and two packets are collided if their 
transmissions overlap in time. Besides, we assume that all DATA packets are of the 
same size for simplicity. We use Scenario 1 for illustration. 
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Scenario 1: 

10 nodes send packets to one AP in a WLAN 
The payload of a DATA packet is 256 bytes 
Fig. 2. Scenario 1 of a wireless LAN 

4.1 802.11 DCF Protocol 



A node with a DATA packet (or simply packet) to transmit first senses channel state. 
If the channel is idle for the time of DIPS (distributed interframe space), the node 
transmits. Otherwise, if the channel is busy during the DIPS, the node backs off, 
i.e., the node defers channel access by a random number of idle slots ranging from 
to CW — 1 within a contention window {CW). When the backoff counter reaches 
zero and expires, the node can access the channel. During the backoff period, if the 
node senses the channel is busy, it freezes the backoff counter and the backoff pro- 
cess is resumed once the channel is idle for a duration of DIPS. To avoid channel 
capture, a node must wait a random backoff time between two consecutive new 
packet transmissions, even if the channel is sensed idle. Once the packet is received 
successfully, the receiver will return an ACK after the duration of SIPS (short inter- 
frame space). SIPS is shorter than an idle slot so that there are no collisions caused 
by DATA packets and ACKs. 

802.1 1 uses the truncated exponential backoff technique to set its CW . In 802.1 lb, 
the initial CW is CWmin = 32. Each time a collision occurs, CW doubles its size, 
up to CWmax = 1024. When the packet is successfully transmitted, CW is reset 
to CWmin- The packet is dropped when it is retransmitted for 6 times and still not 
transmitted successfully. Pig. [3]shows the parameters of 802. 1 lb used in our paper. 

The duration of an ACK is the duration of PHY header plus that of ACK header 
transmitted at basic rate, i.e., '-^"^^Jtf-''^ = 304/^5 ^ 16 idle slots. The duration of 
a DATA packet is the duration of PHY header transmitted at basic rate plus that of 
an MAC header and its upper-layer payload transmitted at data rate. Por example, 
suppose the upper-layer payload is 256 bytes, then the duration of an DATA packet 
is ^ + ^^itS^'^ = 398.5/is ^ 20 idle slots. 



4.2 802.11 Service Curve 



Since equal-sized DATA packets are considered, we measure traffic, service and 
backlog amount in packets in our paper. We measure time duration (e.g. SIPS, 
DIPS, DATA and ACK) in the unit of idle slots and define that one time slot of 
network calculus (t = 1) is equal to L idle slots, where 
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Basic rate 


1 Mbps 


Data rate 


1 1 Mbps 


PHY header 


24 bytes 


ACK header 


14 bytes 


MAC header 


28 bytes 


SIFS 


10 /iS 


DIPS 


50 /iS 


Idle slot 


20 /is 


CWmin 


32 


CWmax 


1024 



Fig. 3. 802.11b parameters 



L = {DIFS + DATA + SIFS + ACK) in idle slots. (41) 

Ideally, an 802.11 node transmits 1 packets per time slot (L idle slots in length). 
Suppose the DATA payload is 256 bytes, L = 3 + 16 + 20 = 39 idle slots. 

Sometimes in the paper, "idle slot" refers to ttte time period which equal to ttte 
length of an idle slot and it may not be idle. To avoid this confusion, we will use 
''idle slot" (italic) to denote that the "idle slot" is indeed idle. 

An 802.1 1 node can be modeled as an ideal server (1 packet transmitted per time 
slot) with the impairment process / which is due to contention with the other nodes 
in a WLAN. In practice, it is difficult to calculate / accurately since / depends on 
the complex interactions of traffic arrival and DCP. In this section, we assume the 
saturated state and use A. Kumar's fixed-point model of 802.1 1 ||30l . This model is 
justified to be very accurate in practice [i3T1l . 

We assume that the system is working at the saturated state, that is, the backlog at 
each node is always nonzero. Por a given node, let r denote its transmission attempt 
probability per idle slot and let r] denote the conditional collision probability when 
it transmits a packet. We assume i] is constant and independent for each transmis- 
sion. Intuitively, this assumption becomes more accurate when the number of nodes 
n increases. In ll30ll . the authors derived two general formulas relating r to t]. The 
first one is 

^ _ ^+V + V^ + ■■■+V^ ^42) 

^0 + vbi + v'^b2 + ... + ri%e 
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This equation can be explained as follows. The numerator is the expected number 
of transmission attempts of a packet. The denominator is the expected total backoff 
duration (in idle slots) of a packet, where 6j is the mean backoff duration after 
the ith collision plus 1 (the 1 refers to the first idle slot of a packet transmission), 
hi 802.11, bi = where < i < 6. A packet suffering 6 consecutive 

collisions will be dropped from its buffer. In our calculations, we do not consider 
packet drops. Since the probability of packet drops is very small, this simplification 
relaxes the backlog bounds very slightly. 

The independence assumption of t] implies that each transmission sees the system 
at steady state. Therefore, each node transmits with the same probability r. This 
yields 

77 = 1 - (1 - r)"-\ (43) 
Combining Eq. (|42|) and Eq. (l43l) . we can solve r and rj. 

We introduce the following terms. The probability of no transmissions at an idle slot 
in the WLAN, denoted by P„(, is (1 — r)". The probability of having at least one 
transmission at an idle slot in the WLAN, denoted by Pt, is 1 — Pnt- The probability 
of a given node starting a successful transmission at an idle slot, denoted by P,, is 

r{l-v)- 

Fig. Hlplots Eq. (l42l) in dashed line and Eq. (l43l) in solid line when n = 10, 20 and 
100. The intersecting points are the solutions to 7] and r. It can be seen from the 
figure that i] increases and r decreases as n increases. Consequently, Pnt increases 
while Pg and Pt decreases as n increases. When we consider the saturated state of 
the system, we actually consider all nodes contending the channel which gives the 
worst-case analysis of the impairment process of a given node and thus conservative 
backlog bounds of it. However, we argument that it is necessary because one applies 
network calculus to deriving the worst-case bounds. 

In order to characterize the impairment process / of the given 802.11 node, it is 
crucial to know its moment generating function. Specifically, we want to calculate 



M7(t) = sup[Ee^^(^'"+*)], (44) 

s>0 

and then we can know its {o-(9), p(6'))-upper constrained characterization (we re- 
place 02 by 6 here to simplify explanations). 

We calculate Eq.(l44l) for a given 802. 1 1 node as follows. Consider the duration of t 
time slots (i.e., tL idle slots) from stos + t. Since we take the sup, we can assume 
that there is always a transmission by the other nodes at the first time slot [s, s + 1], 
which gives the conservative estimation of /(s, s + t) for the given node. We can 
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Fig. 4. The Plots of Eq. (02]) and Eq. ^ 



see that it is a good approximation as the transmissions by the other n — 1 nodes 
happen much frequently than the given node when n is large and also an idle slot is 
much smaller in length than L (in other words, the channel is often busy). We can 
see this point is right for Scenario 1 in the end of this subsection. 

In the following, we consider probabilistic events in the remaining t — 1 time slots. 
There are two cases. Case I: The last transmission is "incomplete". Case II: Other- 
wise to Case I. By "incomplete", we means that the last transmission goes on for k 
idle slots (l < k < L — 1) and get truncated due to the boundary of the last time 
slot. Let Ps = Ps/Pt denote the condition probability of the given node's successful 
transmission on the condition that there is a transmission on the channel. 

We first calculate M/(t) for case I. Suppose there are i complete transmissions and 
one incomplete transmission occupying k idle slots, its probability denoted by pi^k 

is Pt ■ C^(^t-i-i)L-k+i^tPnt~^~^^^~'' ■ here, we use the fact that there are i complete 

transmissions, 1 incomplete transmission and thus (t — i — 1)L — k idle slots in 

the remaining t — 1 time slots. Suppose there are j successful transmissions from 

the given node in the i complete transmissions, its probability denoted by j 

C/ (Ps)-' (1 — PsY'^ ■ When the last incomplete transmission is from the other nodes, 

Miit) is e^^^~^^; otherwise, the transmission is from the given node itself, Mj(t) is 
ge(t-j-fc/L)_ 

Numerating all possible k, i and j, Mj{t) under case I is 
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k=l 1=0 j=0 

= T.Y.P^APse-''/'' + 1 - Ps){Pse-' + 1 - PsTe'' 

k=l i=0 

(45) 



Then we calculate Mj{t) for case II. Let Pi = C}^_^_^-^^_^^PlPnt * ■ Following 
the similar arguments as above, Mj{t) for case II is: 

E Ep.?..e^^*-^'^ = Y^P^iPsC-' + 1 - P,ye''. (46) 
Adding Eq.(l45]) and Eq.(|46l), finally we get Mi{t). 

In general, we do not have the analytical form of Mj(t), so we resort to numerical 
methods to obtain ai{9) and pi (9) (see Algorithm 1 in Appendix B). The algorithm 
is immediately inspired from Definition [6l Then we can use Eq. (I22l) and (|23] ) to 
obtain the node's stochastic service curve. 

We illustrate the above calculations for Scenario 1 in Fig. [2l From Eq. (|42l) and 
dH, T = 0.037 and r] = 0.293. Thus, P^t = 0.680, Pt = 0.320 and = 0.027. 
Again, we can see that the previous assumption that the first slot in [s, s + t] is 
occupied by a transmission from the other n — 1 node is a good approximation, as 
the transmissions by the other n — 1 nodes happen much frequently than the given 
node when n is large (compare Pt — Ps and Ps here) and also an idle slot is much 
smaller in length than L (39 idle slots here; in other words, the channel is often 
busy). 

Fig. |5] shows J's cr{9), (p(6'))-upper constrained characterization when 9 ranges 
from 0.01 to 5.0. 

For example, when 9 = 0.1, we have (T/(0.1) = 0.077 and p/(0.1) = 0.924. And 
(3{t) = {l-rj)-t 

„0.0077 

a(x) = - 6^°-^"^ 

' I _ g0.0924-0.lr/ ^ ' 

for any 1 > r/ > 0.924. 

Finally, we can not apply Proposition[T]to tightening g{x), as p(9) + a{9) > 1 here 
(considering Mf (1) = < e^^^^^"'^^'^). In order to apply this proposition, we must 
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(a) 




Fig. 5. The impairment process I's {a{9), p(0))-upper constrained characterization 

let ri > p{e) + a{e) > 1 which makes /3{t) = (1 - rj)t < 0. 



4.3 Arrival Curves 



In our performance evaluation we use Poisson traffic and we let A be the average 
rate (packets/slot) of it. We have = A by Definition |7l and 



^^eA{s,s+t) ^ j2 ^ ^ c'^^c^' = e^*^'''"^^ (47) 
i=0 ^' 

. Therefore, Poisson traffic is ((Tyi(^^), p^(6')-upper constrained where caI^) = 
andpA = ^. 

We can get the vb arrival curves by Eq. (f2T)) . Note that we can improve the bounding 
function by f{x) = e^^^ for > Pa{0) by Proposition [T] as Poisson process is 
time-independent . 

As for the traffic in reality, we can get Ee^^('*''^+*^ from traffic traces. Then we use 
Algorithm 1 to get the ((Ta(^)) Pa (^) -upper constrained characterization and the vb 
arrival curves. 
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4.4 Stability Condition and Backlog Bounds 



By Proposition m the stability condition of an 802.11 node in a WLAN is 



The stability condition in Scenario 1 is < 0.079 packet/slot or 0.207Mbps by 
the 802. 1 1 parameters in Fig. [3l 

The backlog bounds is calculated immediately by Eq. (l24)) by plugging into traffic 
arrival curves and the 802.11 node's service curve. Note that it is an optimization 
problem depending on 6i and 62- In general, we do not have an analytical solution 
for it. Since the problem dimension is very small, we can apply the method of 
exhaustion to get the optimal value. 

We can also use Theorem |7] to improve on backlog bounds, as in our model we 
consider I under the saturated state which is independent of A. And it still needs to 
optimize the derived bounding function. 



5 Performance Evaluation 

In this section, we compare our backlog bounds derived in Section |4] with ns-2 
simulations in Scenario 1 with Poisson traffic arrivals. The duration of each ns-2 
simulation is 100 seconds which is long enough to let a node transmit thousands of 
packets. And we get the real P{B(t) > x}) over 100 independent simulations. 



As shown in Section 14.41 we can derive stable backlog bounds when = A < 
0.079 packet per slot. Fig.[6]plots the average backlog E[i?(t)] of ns-2 at t = 50s 
and A = 0.077, 0.079 and 0.081 packet/slot. We note that there is a sudden jump 
when A = 0.081, indicating the critical point of stability is indeed around 0.079. 

We use Jiang's model to calculate backlog bounds since the two models are proved 
to be equivalent in our studied case. There are two results of vb arrival curves for 
the traffic arrival process A: the general one (Theorem [5}s Eq.(fT8l)) and the time- 
independent one which improves on the former for the time-independent A (Propo- 
sition [T]). And there are two relations between A and the impairment process / of 
the server: the general one (Theorem [2l) and the independent one which improves 
on backlog bounds for independent A and / (Theorem |7]). 

So we can generate four results of backlog bounds. 



< 1 — o/ 



Ps-L 



(48) 



Put + PfL 
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0.076 0.077 0.078 0.079 0.08 0.081 0.082 0.083 
avg. arrival rateX 

Fig. 6. EB{t) (t = 50s) when A = 0.077, 0.079, 0.081 

• Bound 1: by general vb arrival curves of traffic and general backlog bounds. 

• Bound 2: by time-independent vb arrival curves of traffic and general backlog 
bounds. 

• Bound 3: by general vb arrival curves of traffic and A-I -independent backlog 
bounds. 

• Bound 4: by time-independent vb arrival curves of traffic and A-I -independent 
backlog bounds. 

Obviously, Bound 1 is the loosest and Bound 4 is the tightest among them. 

To illustrate our results, we show the smallest x that makes P{B{t) > x} < p for 
someprobability p, i.e., min{x : P{B{t) > x} <p} where]9 = 0.9, 0.8, 0.1, 0.05 
Obviously, for the same p, smaller x, tighter the bound. Fig. |7] (Fig. [8]) shows our 
results when the 802.1 1 node's traffic arrival rate is A = 0.04 (0.07) packet/slot. 

We make the following observations. First, the backlog bounds improve signifi- 
cantly when we apply time-independent vb arrival curves for traffic or A-I inde- 
pendent case analysis. Second, the bounds of network calculus are much looser for 
higher traffic arrival rate while the real bounds of ns-2 simulations do not relax 
much. Note that A. Kumar's 802. 1 1 model become very accurate near the saturated 
state [.31,1 which is just the case here. The actual reason is: We have the constraint 
of + P/(^2) < c (Eq. (|24l) ). And we need to make 9i and 62 smaller to satisfy 

this constraint for higher traffic arrival rate, which leads to much looser bounding 
functions. Moreover, Theorem[5}s Eq. (fT8l) applies Boole's inequality to the bound 
functions of /, which are loose in general. Here brings the challenge for better 
network calculus models. Finally, we found in trace files that backlog bounds are 
sensitive to the parameters (i.e., 9i, 62, ta and r/) and it is necessary to optimize 



25 





min{a; : P{B{t) > x} < p} 


p 


Bound 1 


Bound 2 


Bound 3 


Bound 4 


ns-2 


0.9 


24 


8 


7 


7 




0.8 


25 


9 


8 


8 




0.7 


25 


10 


8 


8 




0.6 


25 


10 


9 


9 




0.5 


26 


11 


10 


10 




0.4 


27 


12 


10 


10 




0.3 


28 


13 


11 


11 




0.2 


29 


14 


12 


12 


2 


0.1 


31 


17 


14 


14 


2 


0.05 


33 


19 


16 


16 


4 



Fig. 7. ns-2 and network calculus results of backlog bounds under A = 0.04packet/slot (low 
traffic load) 
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Fig. 8. ns-2 and network calculus results of backlog bounds under A = 0.07packet/slot 
(severe traffic load) 

them. 
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6 Related Work 



In this section, we first present a brief overview of the theories and applications of 
stochastic network calculus and then the related works on the performance analysis 
of 802.11. 

The increasing demand on transmitting multimedia and other real time applications 
over the Internet has motivated the study of quality of service guarantees. Towards 
it, deterministic and stochastic network calculus has been recognized by researchers 
as a promising step. 

Essentially, the network calculus is the theory of queueing systems that comes from 
the seminal work by Cruz on the (cr, p) traffic characterization Q []3l and work on 
the service curve characterization of Generalized Processor Sharing (GPS) sched- 
ulers flU [|5l. The theory has been developed by many researchers since then. The 
elegance of network calculus is due to the fundamental convolution formulas (under 
the min-plus algebra) that determine the departure process of a system from its ar- 
rivals and its service curve. The notable strength of the min-plus convolution is the 
ability to concatenate tandem nodes along a network path, and therefore network 
calculus has the ability to characterize the whole network as a single server, which 
is generally intactable by traditional queueing theory OJ. Le Boudec's book cov- 
ers deterministic network calculus and its applications in the Internet [6J. Chang's 
book substantially presented the first approaches to stochastic network calculus 
besides deterministic network calculus fT\. Jiang summarized different types of 
stochastic arrival and service curves in a unified framework and proposed a new 
stochastic network calculus model stemmed from mb (maximal backlog centric) 
arrival curves, although its application conditions have some unsolved controversy. 
Jiang also wrote a book on the theory of stochastic network calculus [8|. Ciucu 
proposed an effective stochastic service curve that can be applied to concatenated 
systems and calculating end-to-end delay and backlog bounds, which exhibits a 
good scaling property of 0{H\ogH) where H is the number of nodes traversed 
by a flow [fTTl . Ciucu also showed that his model can derive quite accurate de- 
lay bounds in M/M/1 and M/D/1 queueing systems by using the martingale tech- 
nique [|T8l . More recently, Fidler proposed a novel solution of the queue system 
using expectations instead of probabilities ||25l , and he also made a comprehensive 
survey on the recent progress of stochastic network calculus [9|. Besides, Jiang 
wrote an overview on this topic from the queueing principle perspective and he 
presented a nice outlook by discussing many open challenges [fTOl . 

Many works have applied network calculus, for example, in measurement-based 
admission control schemes |fT9l , in conformance testing, EOl . in wireless sensor 
networks [i21il . in Aloha systems |I26J, in speeding up network simulations ll22l 
[|23l . in bandwidth estimation ll24l and even in manufacturing blocking systems in 
management science l|27l . 
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Compared with the existing theories of stochastic network calculus, we study the 
effectiveness of backlog bounds in a practical 802.11 WLAN by using the two 
classic stochastic network calculus models: Jiang's and Ciucu's. The latter can be 
extended to concatenated systems while the former still can not at the moment. 
Interestingly, we find that the two models can derive equivalent stochastic service 
curves and backlog bounds in our studied case, which can provide some hints for 
unifying the theories of stochastic network calculus in the future. 

Existing works on the performance of 802. 1 1 focus primarily on the throughput and 
capacity. Bianchi proposed a Markov chain model of 802.1 1 [29 1. A. Kumar et al. 
proposed a probability model of 802. 1 1 [[30ll which simplifies Bianchi's model and 
it is shown to be quite accurate even in the multi-hop case [31J. In our paper, we 
adopt A. Kumar's model to derive the stochastic service curves of an 802.1 1 node. 
There are some works on 802.11 queueing analysis based on traditional queueing 
theory. Zhai et al. assumed Poisson traffic arrival and proposed an M/G/1 queueing 
model of 802.11 (SH. Tickoo proposed a G/G/1 queueing model of 802.11 (SSl 
[[34]l . Bredel and Fidler modeled the 802. 1 1 DCF as a fluid GPS scheduler yielding 
a fair average service rate [25J. And Ciucu analyzed the non-asymptotic throughput 
and delay distribution in multi-hop wireless networks by network calculus approach 
considering Aloha systems [|26l . 

Compared to existing analysis of 802.11, we are the first to analyze the concrete 
802.11 transmissions by Jiang's and Ciucu's models and study their ejfectiveness 
on bounding backlogs. We show that stochastic network calculus is effective theo- 
retically in that the bounds imply stable backlogs as long as the average arrival rate 
is less than the average service rate. However, the bounds are quite loose and we 
show that they can be improved significantly for time-independent arrivals or under 
the independent cases of arrival and service processes. And we note that it is still a 
challenge for a better theory of stochastic network calculus towards tighter bounds 
in practice. Therefore, our work offers a good reference to applying stochastic net- 
work calculus to practical scenarios. 



7 Conclusion and Future Work 

In this paper, we present concrete computations of 802.11 backlog bounds and 
study the bounds effectiveness using stochastic network calculus, from general 
models to detailed calculations. We model a wireless node as a single server with 
impairment service based on two best-known models in stochastic network cal- 
culus: Jiang's |fT6l and Ciucu's IfTTll . And we find that they can derive equivalent 
stochastic service curves and backlog bounds in our studied case. Then we care 
about the effectiveness of network-calculus backlog bounds theoretically. And we 
prove that the network-calculus backlog bounds imply stable backlog as long as 
the average rate of traffic arrival is less than that of service. Next, we consider 
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the effectiveness of network-calculus bounds in practice. We derive the stochastic 
service curve of an 802.11 node from A. Kumar's 802.11 model, which is crucial 
to get backlog bounds. We observe the derived bounds are loose when compared 
with ns-2 simulations. However, the martingale and independent case analysis tech- 
niques can improve the bounds significantly. But still the bounds are not tight. We 
note the reason is due to the looseness in network calculus itself such as Theo- 
rem's Eq.(fT8l). The open questions are: Can we find tighter bounding functions 
under certain conditions? How do we optimize on stochastic arrival/service curves? 
Furthermore, does there exist any unified theory of all network calculus models? 
And these are the future works. 
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Appendix A: Proof of Proposition [T] 

Proof: 

For a fixed t, we construct a stochastic process X(s) = Q^iMi-'^,'t)-rs) (0 < s < t) 
and we have X{s + 1) = X{s)e^^^^^-'-^'*-''^-''\ We will show that if a{s) = 
A{s) — A(s — 1) is independent for each time slot s and r > p{6) + (7{6), then X(s) 
is supermartingale, i.e., E[X(s + 1)|X(0), X{s)] < X{s). 

Because a(s) is time-independent, we have 



E[X(s + 1)|X(0), X{s)] = X{s) ■ E[e^("(*-^)-^)] 

= ■ e"^" ■ Ee^"(*-'\ (49) 

Because A(t) is ((t(6'), p(6'))-upper constrained, we have Ee^'^'^'^^ < e^^^^^'^^^^ for all 
s > 0. When r > p(6') + a{e) and by Eq. (gH), we have 
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E[X{s + l)\X{0),...,X{s)]<X{s). 



(50) 



Thus, is a supermartingale. 

Doob's martingale inequality says that P{suPq<3<(X(s) > A;} < ^^^^ whenX(s) 
is a supermartingale (note: EX(0) = 1 here) for any constant k. Let k = e^, we 
have 

P{ sup [{A{s,t) - r{t - s)] > x} 

0<s<t 

= P{ sup [X{s)] > e^} < e-^ (51) 

0<s<t 



Appendix B: Algorithm 1 (Numerical Calculation of ai{6) and pi{9)) 

Let y{t) = supg>o{| logEe^^'^''''*+*^}. Obviously, y(t) is an increasing function of t 
with y{0) = 0. We define axes t and axes t± (vertical to t) on a plane and we can 
imagine plotting y{t) on it. We define the slope of y{t), s{t) = y{t) — y{t — 1). 

We calculate s{t) for t = 1,2, 3, ... until it converges at some t*, i.e., (1 — e) ■ s{t* — 
1) < s{t*) < (1 + e) ■ s(t* — 1) where e is a small number, e.g. 10"^. 

We draw a straight line l(t) with the slope s{t*) crossing the point (t*,y(t*)^ on 
the axes of t and t_i_. Obviously, the line crosses the point (o, y{t*) — s{t*)t*^ . The 
maximum vertical distance between y{t) and l{t), Vm = maxo<t<t*{y(t) — l{t)}. 
We shift l{t) by and get l{t). Clearly, y{t) < l{t). By Definition^ p{e) = s{t*) 
and a{9) = y{t*) - s{t*)t* + Vm- 
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